source("funcs_and_cons.R")

## production function

pf <- fread(paste0(dir.generated,'matlab/Calibration/',calibration_version,'_KrusellRobotsEach_sol.csv'))

par.kappa_M = rnd2(pf$kappa_M[1])
par.kappa_R = rnd2(pf$kappa_R[1])
par.kappa_C = rnd2(pf$kappa_C[1])
par.kappa_MR = rnd2(pf$nu[1])
par.gamma_E = rnd2(pf$lambda[1])
par.gamma_C = rnd2(1-pf$mu[1])
par.gamma_B_M = rnd2(1-pf$eta_M[1])
par.gamma_B_R = ceiling((1-pf$eta_R[1])*100)/100
par.gamma_B_C = rnd2(1-pf$eta_C[1])
par.gamma_M = rnd2(pf$upsilon[1])
par.A = rnd2(pf$A[1])
par.alpha = rnd2(pf$alpha[1])
par.rho = rnd2(pf$rho[1])
par.sigma = rnd2(pf$sigma[1])
par.nu = rnd2(pf$scale_L[1])
